Bi-quadratic magnetoelectric coupling in underdoped La2Cu04 +a; 
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The recent discovery of relaxor ferroelectricity and magnetoelectric effect in lightly 



5/3 ■ doped cuprate material La2Cu04+ x has provided a number of questions concerning 



its theoretical description. It has been argued using a Ginzburg-Landau free en- 
ergy approach that the magnetoelectric effect can be explained by the presence of 

ergy functional, we study the variety of behavior which can emerge in the electric 
polarization under an external magnetic field. Subsequently, we discuss the role of 

> 

0^ ' Dzyaloshinskii-Moriya interaction in generating this magnetoelectric response. This 

00 ; 

work is particularly relevant for such relaxor systems where the material-dependent 
I/"") ! parameters would be affected by changes in e.g. chemical doping or cooling rate. 
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"jj ! INTRODUCTION 

The parent cuprate materials are Mott insulating antiferromagnets that display a wide 
variety of novel ground states upon doping, including glassy magnetic phases and uncon- 
ventional superconductivity. However, a detailed understanding of how the dopants lead to 
such complex physics is yet to be obtained. Some open issues are related to the structural 
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properties of the dopant atoms such as their preferred locations within the host lattice, the 
physics of e.g. oxygen ions to form clusters, and also the dynamical properties of dopants. 
One interesting example in this context is the oxygen doped material La2Cu04 +x (LCO) 
which, in the undoped case, contains an antiferromagnetic phase with Neel temperature 
TV ~ 325K. Structurally the oxygen dopant ions take spatially non-stoichiometric posi- 
tions unlike e.g. strontium that acts as a substitutional dopant by replacing lanthanum ions 
from the host lattice. Experimentally, a small amount of non-stoichiometric oxygen ions will 
nearly always be present in LCO which, as shown in Refs. 1, 2], can be utilized to study the 
unusual physical properties of the very lightly doped cuprates. 

We have recently discovered that the nearly undoped La2Cu04 +x is a relaxor ferroelec- 
tric at low temperatures T 1 . jl| Further, there exists a weak magnetoelectric effect that is 
anisotropic for different directions of electric polarization and external magnetic field. In a 
recent paper, we were able to reproduce the qualitative features of the magnetoelectric curves 
using a Ginzburg-Landau (_GL) analysis that was governed by the effects of a bi-quadratic 
magnetoelectric coupling. |2| Here, we explore the parameter space of the GL free energy 



re predicted enhancement of 
It is important to analyze 



to study additional effects in the behavior of P(H) and in t 
magnetization and magnetocapacitive effect proposed in Ref. 
these variations in the coefficients of the magnetoelectric coupling particularly for relaxor 
systems, since the coefficients of the GL theory are generally material dependent and may 
be affected by experimental conditions such as e.g. the concentration of dopants, or cooling 
rate. 

MODEL 

The magnetic structure of LCO is that of a two-dimensional antiferromagnet with weak 



feel order. {3- 5| 

The Cu spins are slightly canted out of the Cu02 planes because of a finite Dzyaloshinskii- 
Moriya (DM) interaction existing in the low temperature orthorhombic phase (LTO). 6] On 
application of an external magnetic field, a first order spin flop transition is observed at a 
critical magnetic field H s { ~ 5T. Clear evidence for coupled spin and charge degrees 
of freedom in these systems come from the observation of pronounced discontinuities in 



resistivity and dielectric constant at a magnetic field corresponding to if s f- 
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Below, we follow the theoretical model introduced in Ref. I2I and study the parameter 
space for the magnetoelectric interaction terms that are responsible for spin-charge coupling 
in LCO. We further look at the feedback effect on magnetization at low T due to the presence 
of electric polarization. Below approximately 530K the crystal structure of LCO is LTO with 
space group Cmca (D^). Taking into account the symmetry properties of the Cmca space 
group, the free energy can be expressed as a sum of three contributions 

F = F M + F MP + F P . (1) 

Here, Fu is the magnetic free energy, Fmp is the magnetoelectric contribution, and Fp is the 
polarization free energy. The magnetic free energy that accounts for the crystal structure of 
the LTO phase has been studied previously by e.g. Thio et al. [9| and is given by 

Fm = ~ Ll + \AL\ + \bL\ - CL lMl 

1=1 

+ - H c Mi ~ H ab U] + i J ± LxL 2 . (2) 

Here, the out-of-plane (c direction) [in-plane (a — b plane)] applied magnetic field is repre- 



sented by H c [H a i,]. The coefficients A, B, and C are in general T dependent. [101] The order 
parameter Mj = (Sa% + Sbi)/2 is the ferromagnetic moment per spin with SAi,Sp,i being 
the sub-lattice spins in the i th plane, and = (Sai — Sbi)/2 is the antiferromagnetic order 
parameter (Lj||a). The spins are slightly canted due to the DM interaction term —CMiLi, 
which causes them to lie in the a — c plane of the magnetic unit cell. The coupling between 
the different planes is included by the J± term. 

The presence of an inversion symmetry in the space group of the crystal forbids any 



linear magnetoelectric effect and the physics is dominated by non-linear coupling terms. 
We can focus on the largest non-linear terms by further noting that the experimentally 
observed polarization response is symmetric under inversion of the external magnetic field 
(i.e P{H) = P(—H)). This implies that the dominant couplings are of even order in the 
magnetic order parameter. Hence, the following terms contribute to the magnetoelectric 
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coupling 



MP 



^L'i + ^M'f + ^M^Pl (3) 



where the components for P run over a = (a, b, c) in the magnetic unit cell. 
The polarization free energy is given by 

F p = Dlf P ' + I P «) " EP ( 4 ) 

a 

Here, Xea is the electric susceptibility for the a component of the polarization, and E denotes 
the applied electric field. The solutions that determine P(H) are obtained by minimizing 
F with respect to the electric polarization and magnetic order parameters. In the case of 
LCO studied experimentally, XJv ~ 320K, which is much higher than the T at which the 
ferroelectric order sets in (Tp ~ 4.5K).jl(] Therefore, we evaluate Fm for the high T phase 
with P = 0, providing the following set of equations 

M^xoiHc + CL,), (5) 
XznL-i + ALi + BL\ + l - J ± L 2 = CM, + H ab , (6) 
X ^L 2 + AL\ + BL\ + l -J L L Y = CM 2 + H ab , (7) 

2 

[X« + E(^ L * 2 + + 7 3a M J L l )]P Q = -PPI (8) 

i=l 

The experimental magnetization curves at low temperatures (<~ 30K) have a glassy 
contribution. Q These features cannot be obtained from the above equations and to include 
them to lowest order, we take the experimental magnetization values at T = 5K as input to 
the model. 



RESULTS 

In terms of the following rescaled quantities l + = xoC(L\ + L 2 )/2, /_ = xoC(Li — L 2 )/2, 
M — (Mi + M 2 )/2, i la = 2 7la (xoC)- 2 , i 2a = 2 72a , j' 3a = 2 l3a ( Xo C)-\ the polarization 
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FIG. 1: (Color online) Theoretically calculated P c (H c )/P c (0) for different values of 
Parameters not mentioned in the plots have been fixed at s c = 0.1, s a = 0.074, g a = 9b = 
q c = —6.58, and q a = Qb = that explain the experimental curves at T = 5K. 



dependence on the applied magnetic field can be expressed as 



+ (1 - 9a ~ q a )l 2 + (H c ) + q a M{H c )l + {H c ))} l /\ 



q c and g c . 
9c = -0.2, 



(9) 



(10) 



where s a = \Jl (0) /(Xed + ^J- (°))> 9a = ikJ X, and q a = 7 ^/A Q with \ a = y ia +l' 2a +lL- 
In general, all three parameters s a , g a , and q a are T dependent. The T dependence of s a 
results primarily from its relation to the electric susceptibility. 

In the case of an in-plane magnetic field H a b, P a (H a b) depends on a single fitting param- 
eter s a that only controls the magnitude of the polarization ratio, whereas the shape of the 
theoretical curves are governed by the magnetic order parameter of the system. Q 

In the case of an out-of-plane magnetic field H c , the measured P C (H C ) increases with 
magnetic field for positive values of q c and exhibits a pronounced hump at the spin-flop 
transition at H s f as seen from Fig. [TJ It can be noted that the qualitative features of the 
observed experimental result jl|, is therefore reproduced for an attractive DM induced 
magnetoelectric coupling q a P^MiLi as the polarization would be suppressed with increasing 
field for repulsive DM induced term. Comparing the magnitude of the polarization ratio 
at q c = in Fig. [T] with the value that qualitatively matches with the experiments at 
q c = —6.58 we find that there is a dominant contribution from DM induced magnetoelectric 
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coupling to the polarization enhancement. On the other hand, we can also see from Fig. [TJ 
that the coefficient of bi-quadratic magnetoelectric coupling between the magnetization and 
the electric polarization given by g c does not have any significant influence on the P C (H C ). 

For polarization values measured in the Cu02 plane, the effect of the parameters on the 
P a (H c ) curve is shown in Fig. [2] As can be seen from the figure, an increase of the attractive 
q a term leads to the suppression of the electric polarization in contrast to the behavior of 
P C (H C ) in the presence of finite q c . Though the experimental curve of Ref. [l is reproduced 
by q a — 0, the qualitative features of the curve are maintained for a finite q a <= 0. It is also 
interesting to note the effect of the magnetization coupling term g a on P a (H c ). Unlike the 
case for out-of-plane polarization discussed above, this term strongly influences the form of 
the curve and we find significant deviations from the qualitative shape of the experimental 
curve for different values of g a . It is therefore quite remarkable that the value of g a that 
reproduces the experimental curve corresponds to g a = g c = —0.2 since we do not have 
any a priori reason for it to do so. Also note that the deviation in the polarization are 
much stronger above the spin flop transition. This is to be expected since the magnetization 
undergoes a significant enhancement above the transition. 

We have observed experimentally that the magnetization shows a small upturn below the 



temperatures where the ferroelectric order sets in. [12J This effect is in addition to the typical 
upturn in magnetization near the spin glass freezing temperature. {13I. 1^ | The inclusion of 
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FIG. 2: (Color online) Theoretically calculated P a (H c )/P a (0) for different values of 
Parameters not mentioned in the plots have been fixed at s c = 0.1, s a = 0.074, g a = 9b = 
q c = —6.58, and q a = qb = that explain the experimental curves at T = 5K. 



q a and g a . 
9c = -0.2, 
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such a feedback effect leads to the following expression for the magnetization, 



X qH c + [i - xo E a jLPjmMHo) 

, [1 - XoE a lLP*(H ab )]l + (H ab ) 

1 + 



(11) 



(12) 



Note that in this expression the relative sign of the coefficients can be determined from the 
relations j'^/j^a, = q a /g a > 0. Therefore for values of q a and g a that are opposite in sign 
we would find a suppression in magnetization and an enhancement otherwise thus providing 
another experimental test to identify the relative sign of the GL coefficients. Similarly, 
since the magnetocapacitive effect is proportional to the polarization enhancement due to 
an external magnetic field the size of the jump would be indicative of the magnitude of the 
DM induced magnetoelectric coupling. Q 

In summary, we have studied the parameter dependence of the recently discovered magne- 
toelectric effect in extremely underdoped L^CuC^+a; modeled by a GL free energy including 
a bi-quadratic magnetoelectric coupling term. Changes in this term can lead to interest- 
ing behavior of the field dependence of the polarization, particularly above the spin flop 
transition. 
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